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' Motivated by experiments on the layered compounds k-(BEDT-TTF)2X, CS2CUCI4, and very 

£Sj ' recently Na^OoC^ • 1/H2O, we present a weak-coupling functional renormalization-group analysis of 

the Hubbard model on the anisotropic triangular lattice. As the model interpolates between the 
nearest-neighbor square lattice and decoupled chains via the isotropic triangular lattice, it permits 
the study of competition between antiferromagnetic and BCS Cooper instabilities. We begin by 
reproducing known results for decoupled chains, and for the square lattice with only nearest-neighbor 
£ — ■ hopping amplitude t\. We examine both repulsive and attractive Hubbard interactions. The role 

of formally irrelevant contributions to the one-loop renormalization-group flows is also studied, and 
these subleading contributions are shown to be important in some instances. We then observe 
that crossover to a BCS-dominated regime can occur even at half-filling when antiferromagnetism is 
frustrated through the introduction of a next-nearest-neighbor hopping amplitude ti along one of the 
two diagonal directions. Stripes are not expected to occur and time-reversal breaking d x i_ y i ±id xy 
superconducting order does not arise spontaneously; instead pure d x 2_ y 2 order is favored. At the 
isotropic triangular point (ti — tz) we find the possibility of re-entrant antiferromagnetic long-range 
order. 
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. PACS numbers: 74.20.Mn, 74.25.Dw, 74.70.Kn, 71.10.Fd 
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I. INTRODUCTION 



The behavior of strongly correlated electrons moving in reduced spatial dimensions continues to yield surprising new 
physics. For example, intriguing experiments Q, 0, on the k-(BEDT-TTF)2X family of layered organic molecular 
crystals evoke similar findings in the field of high-temperature cuprate superconductivity. Like the high-T c cuprates, 
the layered organic materials exhibit a wide variety of electronic properties. In particular, the phase diagram is rather 
similar to that of the cupratesQ and there is some evidence for unconventional pairing with nodes in the gap from 
NMR relaxation rate[j||(| 0, specific heat EL penetration depth|9l ITo l [ill II2I fTa | . STM spectroscopy [3, mm- wave 
transmission 0| (see however Refs. [l6[ |17[) and thermal conductivity |l8l Il9| measurements. Other experiments, 
however, suggest s-wave pairing Competition between antiferromagnetic and superconducting 

instabilities, seen in the cuprates, also seems to occur in the k-(BEDT-TTF) 2 X compounds [2(j. In contrast to the 
• square Cu02 lattice of the high-temperature superconductors, the organic molecules pair up into dimers, and these 
dimers form a triangular lattice. 

Two other quasi two-dimensional materials with triangular lattices have been the subject of recent attention: 
The antiferromagnetic insulator CS2CUCI4 compound |27j and the cobalt-based superconductor Na^CoO^ • yH 2 that 
may be an analog of the cuprate high-temperature superconductors |28| . Neutron scattering experiments suggest 
the existence of deconfined spinon (spin-1/2) excitations in the CS2CUCI4 once antiferromagnetic order has been 
eliminated by heating the sample to the relatively low temperature of approximately 0.6K, or upon application of a 
field parallel to the plane [2^, Geometric frustration of the spin-spin interactions is likely responsible for the observed 
spin-liquid behavior. As the cobalt atoms in the Na. T CoQ 2 • 2/H2O material also form a triangular latticej^, and as 
they have further been argued to carry spin-1/2 moments |28j. we tentatively group this system into the same category 
as k-(BEDT-TTF) 2 X and Cs 2 CuCl 4 . 

Clearly theoretical investigations of strongly correlated electrons on triangular lattices are of great interest. Initial 
studies of strongly correlated systems often start with a minimal Hubbard model, leaving extensions such as the 
inclusion of long-range Coulomb interactions for later more detailed work. In fact McKenzie has proposed^ that a 
Hubbard model on the anisotropic triangular lattice serves as a minimal model of the conducting layers of k-(BEDT- 
TTF)2X. It represents a simplification of a model introduced earlier by Kino and Fukuvama|3dj . Two distinct hopping 
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FIG. 1: Anisotropic triangular lattice with two hopping amplitudes ti and t^- The limit ti = 0.0 corresponds to the usual 
nearest-neighbor square lattice and t\ — 0.0 corresponds to decoupled chains. 

matrix elements are introduced and the Hamiltonian is defined by: 

H = -t 1 ^{c\ a c- ia +H.c.)-t 2 (4 a ^ + H - c -)+ tf$>iT"ii -^m, (1) 

<ij> <<ij>> i i 

where < ij > denotes nearest-neighbor pairs of sites on the square lattice and << ij >> denotes next- nearest- 
neighbor pairs along one of the two diagonal directions of the square lattice as shown in Fig. ^ Quantum chemistry 
calculations suggest that, unlike the cuprate materials, in the case of the organic k-(BEDT-TTF)2X compounds the 
Hubbard interaction U r* t. Thus a weak-coupling renormalization-group (RG) approach such as we adopt here may 
be expected to be reasonably accurate for the organic materials. 

The model is also interesting in its own right as it interpolates between the square lattice and decoupled chains. At 
half-filling, the non-interacting Fermi surface is perfectly nested in these two extreme limits. As nesting is imperfect in 
between the limiting cases, several phase transitions can be expected. The square lattice, which has been the subject 
of many studies, corresponds to the special case of zero next-nearest-neighbor hopping, t 2 — 0. When the repulsive 
interaction is turned on, nesting induces a spin density wave instability. In the opposite limit, t\ — 0, the chains are 
completely decoupled. These isolated chains of course have no spin order and are described by the exact Bethe ansatz 
solution of Lieb and Wu|3l|. We pay particular attention to the intermediate region of t\ ^ and t 2 ^ and study 
it via a weak-coupling renormalization-group analysis. The special isotropic triangular lattice point corresponds to 
t\ = t 2 . Values for the hopping matrix elements obtained from experiments and from quantum chemistry calculations 
for the conducting layer of k-(BEDT-TTF) 2 X suggest t\ > t 2 , that is, somewhere intermediate between the square 
and the isotropic triangular limits. The lattice anisotropy can be altered by uniaxial stress applied along the principal 
axes of the quasi- two-dimensional organic compound j32|. Fermi surfaces of non- interacting electrons for different 
ratios of the hopping matrix elements are shown in Fig. 

In the next section we briefly introduce the RG method we employ, a method first implemented by Zanchi and 
Schulz 33] for the case of the square lattice. We then present results of our calculations at different values of the 
anisotropy: decoupled chains (studied as a test case to check the reliability of the calculation), square lattice, and 
finally the anisotropic region intermediate between the square lattice and the isotropic triangular lattice. We discuss 
the ordering tendencies and work out the implied phase diagram as a function of anisotropy parameter t 2 /{t\ + 1 2 ) 
which ranges from (square lattice) to 1 (decoupled chains) . We also make comparison to results obtained via other 
methods in the strong-coupling limit of large on-site repulsion. 

II. RENORMALIZATION-GROUP CALCULATION 

We follow the weak-coupling renormalization-group analysis implemented by Zanchi and Schulz [33^ for interacting 
fermions on a two-dimensional lattice. Like some previous work 34], the approach generalizes Shankar's renormaliza- 
tion group theory|3!| to Fermi surfaces of arbitrary shape. More significantly, in principle the only approximation 
that is made in the approach of Zanchi and Schulz is an expansion in powers of the interaction strength, the on-site 
Coulomb interaction U. Subleading terms generated during the RG transformations, which are dropped as irrelevant 
in the simplest versions of the RG, are instead kept in this formulation. Specifically, the formally irrelevant, non- 
logarithmic, terms that appear in the six-point function during the process of mode elimination do in fact contribute 
to the RG flows. Thus while the simplest weak-coupling RG analyses makes a double expansion in both the interac- 
tion strength and in the relevance of the terms retained in the renormalization flows, in the approach of Zanchi and 
Schulz there is only a single expansion in the interaction strength. (In practice some additional approximations are 
made for computational convenience, as detailed below. These simplifications are not expected to alter the results 
significantly.) As we show below, in some cases this more accurate treatment leads to substantial differences in the 
RG flows. 
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FIG. 2: Fermi surface of non-interacting electrons for different ratios of the two hopping amplitudes. The number at the top 
of each graph is the anisotropy ratio £2/(^1 + £2) which ranges from (square lattice) to 1 (decoupled chains). The chemical 
potential fj, is varied to ensure that the system remains half-filled. 

Elimination of high energy modes is carried out iteratively, in infinitesimal steps, and as a result the energy cutoff 
A around the Fermi surface shrinks, see Fig. [31 The initial energy cutoff is taken to be the full band width A , and 
it is reduced via continuous mode elimination to A = Aq e~ l where £ > 0. For each infinitesimal step £ — » £ + d£, the 



Fermi Surface 




FIG. 3: Mode elimination is carried out in infinitesimal steps. The figure shows the low energy modes (denoted L) which are 
inside a shell of thickness A around the Fermi surface, and the high energy modes (H), which have already been integrated 
out and which are now outside the shell. At each step the on-shell modes precisely at scaling parameter £ (dashed lines) are 
integrated out. 



fermion degrees of freedom are broken down into high and low energy modes as 

¥ a (K) = 0(A - |e k |) * a , L {K) + 9(|e k | - A) * a%H {K) , (2) 

where K = (cj, k) is the usual 2+1-dimensional frequency-momentum vector. The effective action, after dropping a 
constant contribution fin to the free energy, has the form 

SA(t+<U) = S&y) + 5S(£). (3) 

At non-zero I > the effective action contains contributions at all orders in the initial interaction strength. But 
because mode elimination is done in infinitesimal steps, only terms linear in d£ contribute to 5S(£). These terms 
correspond to diagrams with one internal line (either a loop or a tree diagram). RG flow equations for vertices with 
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any number of legs, T2 n (£ ), can then be found. These are functional equations since the T's are functions of momenta 
and frequencies. 

To make progress we must make an approximation. We carry out the weak-coupling expansion by truncating the 
RG equations at the one-loop level. Renormalization of the effective interaction L^(ki, k2, k3, IC4), corresponding to 
the four-point function (1^4), then occurs at order U 2 . Contributions from the six-point functions Tq must also be 
included at this order. Higher n-point functions may be neglected as these only contribute at higher-order in the 
interaction strength U. It is important to notice that the RG flow equations generated this way are non-local in 
scaling parameter £. The RG equations for couplings [7g(ki, k2, k3, k^) at step I involve the values of couplings at 
previous steps £ pp and £ p h [the subscript denotes particle-particle (pp) and particle- hole (ph) channels] : 

^—K^*) <4 » 

^ = -l»(^f) (5) 

with q pp = ki + k2 and q p h = ki — k4. At step £ contributions from six-point functions are obtained by contracting 
two of the legs at on-shell momentum k. Of course the six-point functions were generated from four-point functions 
during previous steps. Momentum k4 is determined uniquely by momentum conservation to be k4 = ki + k2 — k3, 
so we drop it in the following. 

For an initial, bare, four-fermion interaction Uq which is independent of spin, following Zanchi and Schulz it 
is possible 33] to write all the renormalized two-particle interactions in terms of only one function Ug (ki, k2, ^3). 
Couplings in the charge and spin sectors can then be obtained from this function through the relations: 

U C = \{2-X)U, U - = ~^U. (6) 

where X is a permutation operator defined by its action: XU (ki, k2, ^3) = U(k.2, ki, k.3). The charge density (CDW) 
and spin density (AF) couplings are then given by: 

Vf DW (9 1 ,9 2 ) = 4[/ rf (k 1; k 2 ,ki) 
V e AF (6 1 ,d 2 ) = 4[/ ff ,(k 1 ,k 2 ,ki) . (7) 

Here kj is related to kj by kj — kj = Q where Q is a nesting vector [Q = (±tt, ±tt) for the fully nested square lattice]. 
Also 9j is the angle that wavevector kj makes with the x-axis. The forward scattering amplitude is given by 

F e (e u 9 2 )^U i (k 1 ,k 2 M) (8) 

and only involves two momenta (ki and k2) because the momentum transfer during scattering is very small. Likewise 
the BCS interaction 

V i BCS (9 1 ,9 2 ) = U e (k 1 ,-k 1 ,k 2 ) (9) 

also is described by just two momenta as it represents the scattering of a Cooper pair of electrons of opposing momenta 
ki and — ki into a pair of electrons of opposing momenta k2 and — k2. 

In order to integrate the flow equations forward in the scaling parameter £ we first discretize the Fermi surface, 
dividing it up into patches as depicted in Fig. Replacing the continuous surface with discrete patches should be 
adequate for the imperfectly nested Fermi surfaces we focus on here [3^. After discretization of the Fermi surface, 
the angles 9j = 2iTj/M where j = 0, . . . ,M — 1. Interactions Ue(ii,i 2 ,is) are thus labeled by three discrete patch 
indices. A further approximation is implied by this procedure, as the dependence of the effective interaction on the 
radial component of momentum is neglected and the shape of the Fermi surface is not renormalized. The justification 
is the following: though the shape of the Fermi surface change at the one-loop level, the feedback of this change on 
the one-loop RG flows for the couplings Ut(i\, 12,13) constitutes a higher-order effect. The dependence of U on the 
radial components of the three momenta is irrelevant |35l 13^] . This is similar to the one-dimensional case, where the 
marginal interactions are labeled according to the indices i = L, R (left or right moving) of the electrons that are 
interacting. There is strong dependence on the direction of k, but the dependence on the absolute value |k| of the 
momentum is irrelevant. Therefore the interactions may be parameterized simply by their projection onto the two 
Fermi points. In two dimensions the interactions are likewise parameterized by the patch indices. 

In this work we only study flows at zero temperature. The integral over Matsubara frequencies, which arises 
in the one-loop diagram, can be performed analytically as the dependence of the couplings on the frequency lj is 



12 



FIG. 4: Discretization of the Fermi surface into M = 16 patches. Each patch corresponds to an angular section of 2n/M. The 
special case of the perfectly-nested Fermi surface corresponding to the nearest-neighbor tight-binding model on a square lattice 
and at half-filling is shown. 

irrelevant [3 7\. We set the initial bare coupling to be \Uq\ = 1 and, unless otherwise stated, also set t\ +t% = 1. The 
full band width is Ao = 8ti +4i2- We usually divide the Fermi surface into M = 16 patches. For the special case of the 
isotropic triangular lattice we instead use a finer mesh of patches, M = 24, to permit an examination of higher-wave 
channels. Our algorithm makes no assumptions about the symmetries of the Fermi surface; this means that we must 
follow the flow of all M 3 couplings U(i\, 12, £3). We do impose the requirement that the three indices are such that 
all four particles lie on the Fermi surface. The RG flow for these couplings are then described by coupled non-local 
integral-differential equations. These equations are numerically integrated forward in the scaling parameter I. The 
increment in the scaling parameter is set to be dt = 0.1 for the results shown here. Calculations using smaller values 
of dl yield nearly the same results. 

An equivalent version of RG method for two-dimensional interacting fermions has been developed by Salmhofer 38, 
l39l |. In this formulation, the RG flow equations are local in the scaling parameter £, but this gain comes at the cost of 
expanding the effective action in Wick-ordered monomials, resulting in RG flow equations with one extra integration 
over momentum. This formulation has been used to study the two-dimensional Hubbard model on a square lattice 
with nearest-neighbor and next-nearest-neighbor hopping amplitudes [4(1 l4ll . 

Given a set of RG flows, we must then interpret the various orderin g te ndencies. One way to do this is by calculating 
susceptibilities towards order, as carried out for instance in Refs. l4Ct l4ll Another approach is to bosonize the fermion 
degrees of freedom, and then determine the ground state of the bosonized effective Hamiltonian semiclassically by 
replacing each boson field with a c-number expectation value. The latter method was adopted by Lin, Balents, and 
Fisher in their treatment of the two-leg ladder system [43^. Klein ordering factors must be treated carefully 0] and 
the resulting weak-coupling RG / bosonization prediction was shown to agree well with the results of essentially exact 
DMRG calculations 45j . We leave the extension of such an analysis to the full two-dimensional problem 34] for future 
work, and make the observation here that in most instances it suffices to simply follow, during the course of the 
RG flow, the most rapidly diverging interaction channel. For instance, the effective BCS interaction V e BCS (ii, 12), 
as defined by Eq. is a symmetric M x M matrix in the patch indices. The various BCS channels are obtained 
upon diagonalizing the matrix. The eigenvector <f> with the largest attractive eigenvalue then represents the dominant 
BCS channel. We also calculate the eigenvectors and eigenvalues of the effective spin coupling V e AF (ii, 12) and charge 
density wave coupling Vf DW 12) to determine the dominant AF and CDW channels. In the following we plot largest 
eigenvalues of the interaction matrices as a function of the scaling parameter, as well as the dominant eigenvectors as 
a function of the patch index, to gain insight into the ordering tendencies. As shown in the next section this way of 
intepreting the RG flows yields the correct physics in the limiting cases of one-dimensional decoupled chains as well 
as the completely nested square lattice. 

III. RESULTS 

We now turn to the results of our RG calculation. We first check the method in the special limiting cases of 
decoupled chains and the pure square lattice. As we reproduce known results in these limits, we then turn to the 
more general problem of the anisotropic triangular lattice. 
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A. Decoupled chains (t\ — 0) 



As a first check, we apply the weak-coupling analysis to the case t\ = and = 1 which corresponds to completely 
decoupled chains. At half-filling, particle-hole symmetry requires /i = and the nesting wavevector is Q = it. 
As quantum fluctuations always suffice to prevent continuous symmetries from breaking in one spatial dimension, 
antiferromagnetic and superconducting order are not possible. Instead the possible phases are classified in terms of 
whether or not charge and/or spin excitations are gapped. Furthermore the three types of marginal interactions are 
often denoted (see Ref. \4(t) spin current (A s ), charge current (A c ), and Umklapp (A u ), where the latter two carry only 
charge and no spin. In terms of our notation we may identify 



A c = 4U C (R,L,R), 
A s = 4U„(R,L,R), 
X u = U(R,R,L) . 



(10) 



As shown in Fig. for repulsive initial interaction (Uq > 0), the spin couplings decrease towards zero, whereas 
the Umklapp and charge couplings diverge in the low-energy limit. This is as expected from the exact Bethe ansatz 
solution |3l| since the system has gapless spin excitations while the charge sector is gapped. On the other hand, for 
attractive initial interaction (Uq < 0) the spin couplings diverge while Umklapp and charge couplings tend towards 
zero. In this case there is a gap in the spin sector and gapless charge excitations. The solid lines in Fig. [^correspond to 
a direct analytical solution of the simple one- loop RG equations for the one-dimensional Hubbard model at half- filling: 



dX s 
~M 
dXc 
dt 
dX u 

dl 



4* 

TT 

-X 2 
—X C X U 

TT 



(11) 



Our numerical solution of the Zanchi-Schulz RG equations agrees quantitatively with the standard one-loop results. 
An exact fit is not expected, because the the Zanchi-Schulz equations also include the renormalization of the charge 
and spin speeds as well as sub-leading non-logarithmic corrections. 




FIG. 5: RG flows of the spin current (A s ), charge current (A c ) and Umklapp (X u ) interactions, in the limit of completely 
decoupled chains (tj = 0). For repulsive (attractive) bare interaction Uo > (Uo < 0), a charge (spin) gap develops and there 
are gapless spin (charge) excitations. Solid lines are the analytical solution of the one-loop RG equations for the one-dimensional 
Hubbard model at half-filling. 
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B. Square lattice (t2 = 0) 

Attractive interactions Uq < induce strong BCS instabilities in the case £2 — of a pure square lattice, in accord 
with expectations. The eigenvector of the dominant attractive BCS channel is plotted in Fig. for the case of 
half-filling, /1 = 0. As expected, the BCS pairing is in the s-wave channel when Uq < 0. From the outset at I = 
the BCS sector dominates all other channels. As the RG flows progress, the BCS channel diverges and thus remains 
the dominant coupling. The same qualitative behavior persists as the system is doped away from half-filling. Fig. [7| 
depicts the RG flows in the dominant AF, BCS and CDW channels both for the half-filled case fi = and away from 
half- filling (/i = 1). Perfect nesting at half- filling drives V CDW to diverge as strongly as V BCS . Away from half- filling, 
there is no perfect nesting so both V AF and V CDW grow at a much smaller rate. In Fig. |S| we compare our results 
to RG flows which include only the leading logarithmic contributions. As the BCS coupling is not driven here by 
the AF fluctuations, there is little coupling between the two channels. Therefore, one-loop RG flows which include 
only the leading logarithmic contributions (diamonds) do not deviate significantly from the more accurate approach 
of Zanchi and Schulz (circles) which includes all the subleading non-logarithmic terms generated at one-loop. 
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FIG. 6: Pairing symmetry of the dominant BCS channel for the case of attractive initial interaction (Uo < 0) on a square 
lattice (£1 — 1 and £2 = 0) and at half-filling (fj, = 0), as revealed by plotting the eigenvector of the BCS matrix with the largest 
attractive eigenvalue. The scaling parameter £ — 2.5. The Fermi surface is divided in M — 16 patches. As the patch index 
goes from n = to n — 15 around the Fermi surface, the angle 9 goes from to 2n. 
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FIG. 7: RG flows of the dominant AF, BCS and CDW couplings for the square lattice with an attractive Hubbard interaction 
(Uo < 0). Plotted are the largest eigenvalues of the corresponding interaction matrices. Solid lines correspond to the half- filled 
case (/j, = 0.0). Nesting occurs at half-filling and V CDW diver ges as strongly as V BCS . Dashed lines correspond to /i = 1.0. In 
this case, nesting is destroyed causing V AF and y CDW to increase at a much smaller rate. 
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scaling parameter / 

FIG. 8: Comparison of results obtained from the approach of Zanchi and Schulz (circles) with those obtained when only 
logarithmic contributions are included (diamonds), for the case of an attractive Hubbard interaction (Uq < 0) and (J, = 1. 
There is no significant discrepancy between the results. Flow in the dominant BCS coupling is shown by the solid line; the AF 
coupling is indicated by the dashed line. 



The opposite limit of a square lattice with repulsive initial interaction Uq > has been extensively studied[33l l40|. 
At half-filling, the Fermi surface is perfectly nested and a strong SDW instability develops. Fig. [5| shows, however, 
that the largest BCS channel, though sub- leading in comparison to the SDW channel, already exhibits d x 2_ y 2 pairing 
symmetry. 
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FIG. 9: Pairing symmetry of the dominant BCS channel for the case of repulsive initial interaction (Uq > 0) on the square 
lattice (ti = 1 and = 0) at half-filling ([i = 0). The different curves correspond to different values of the scaling parameter 
(i — 0.5, 1 and 3). The Fermi surface is divided in M — 16 patches. As the patch index increases from n — to n — 15 around 
the Fermi surface, the angle 8 increases from to 2-rt. 

For the case of a repulsive Hubbard interaction, Uq > 0, we find in contrast to the attractive situation that the 
formally irrelevant terms play an important role. As the initial BCS couplings are all repulsive, Cooper pairing can 
only happen via coupling to the AF channels or via the non-logarithmic corrections to scaling coming from the formally 
irrelevant terms in the six-point functions. We again compare RG flows which include only the leading logarithmic 
corrections (diamonds) against those in which all subdominant contributions at one-loop order are included (circles) 
in Fig. EH for the case fj, = 10~ 4 , that is, slightly away from half-filling. Though qualitatively similar, there is 
considerable quantitative difference. At this small doping, AF tendencies dominate in both cases. Upon further 
increasing the chemical potential, as mentioned above there is a crossover into the d x i_ y -2 BCS regime. The crossover 
occurs much sooner when the subleading terms are included. Finally, at large doping the nesting of the Fermi surface 
is completely eliminated, and neither the AF nor the BCS channels show any strong divergences. 
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FIG. 10: Comparison of results obtained at one-loop level for the case of a repulsive Hubbard interaction and (i — 10 -4 . RG 
flows which include the effects of formally irrelevant terms (circles) are compared with RG flows for which only logarithmic 
contributions are retained (diamonds). Flow in the dominant BCS channel is depicted as a solid line and dominant AF sector 
is indicated by a dashed line. 

C. Square Lattice With Next-Nearest-Neighbor Hopping (t' ^ 0) 

Next we turn to the square lattice with added next-nearest-neighbor hopping amplitude t' along each of the two 
diagonal directions. Weak-coupling RG studies of this Hubbard model have been carried out previously by a number 
of groups^, HE E3 using the formulation of Refs. |H, |3^ We have checked our calculation against these published 
results and find good agreement. The dispersion relation in this case is given by 



ek = — 2ti(cosfc x + cos k y ) — 4t' cosk x cosk y 



(12) 



and Fig. II II show the Fermi surface for the case t\ — —1, t' — 0.05, and fj, = At'. The Fermi surface is centered around 
the r point (±7r, ±7r) and van Hove singularities lie at the Fermi energy 47] . 



71 



-71 




FIG. 11: Fermi surface of the square lattice with nearest-neighbor hopping ti = — 1 and next-nearest-neighbor hopping t' = 0.05. 
The chemical potential /i = 4t' . 

This Fermi surface exhibits the d x 2_ y 2 superconducting instability. It is interesting to go a bit further and address 

the question of whether or not there is spontaneous time-reversal (T) symmetry breaking with the appearance of an 
additional imaginary id xy component to the superconducting order parameter. To answer this question we study the 
implications of RG flows which yield comparable attraction in two channels: one term with d x i_ y i symmetry and a 
second with d xy symmetry. A simple calculation of energetics then suffices to show that the two order parameters 
will phase-align as d x 2_ y 2 + id xy . The standard BCS equation yields a condensation energy of 



AE = E sc -E N = 2J2 

k>k F 



24 



2\A 



(13) 
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FIG. 12: Dependence of the integral I(u) on the real part of it. The maximum of I(u) occurs when Re{u} = 0, that is, when 
u is purely imaginary. 



For couplings V with comparable d x 2_ y 2 and d xy components, the ansatz to maximize the condensation energy 
should be chosen to be A k = 2 _ 2 (k) + u Ad xy (k) , with u encoding information about the relative phase of the two 
components. Substituting this ansatz into Eq. 1131 we may then determine the phase that maximizes the condensation 
energy AE. Fig. 1121 shows the dependence of the sum 



I(u) 



2e k 



k>k F 



A k 



(14) 



on the real part of u. This term is maximized when u is purely imaginary (Re{u} = 0), hence the T-breaking 
pairing symmetry d x 2_ y 2 + id xy is the energetically favored. Physically this is reasonable, as this choice of the phase 
guarantees that a gap forms everywhere along the Fermi surface, lowering the ground-state energy. 

Returning to the square lattice, we find that upon integrating the RG equations for the case ti = — 1, t' = 0.05 and 
/i = it', the dominant attractive BCS channel has d x 2_ y 2 symmetry as expected; see Fig. I13f A). A channel with d xy 
symmetry also appears but it is repulsive in sign; see Fig. H3f B). The strengths of each channel are plotted in Fig. ^] 
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FIG. 13: Dominant BCS channels for the square lattice with ti = — 1, t = 0.05 and /i = 4t , at scaling parameter i = 5. The 
d x 2_ y 2 channel (A) is attractive while the d xy channel (B) is repulsive. The absolute values of the eigenvalues are plotted in 
Fig. 1141 as functions of the scaling parameter t. 

as a function of the scaling parameter I. Since the d xy channel is repulsive, no d xy order will arise, and this may be 
taken as evidence against the formation of spontaneous time-reversal symmetry breaking of the d x 2_ y 2 + id xy type. 
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FIG. 14: Flow of the eigenvalues for the d x 2_ y 2 and d xy channels shown in Fig. ED The magnitude of the eigenvalues is plotted. 

D. Towards the Triangular Lattice (t2 7^ 0, t<z < t\) 

Introducing non-zero ti along just one of the two diagonals, as shown in Fig. ^ offers a different way of breaking- up 
perfect nesting and enhancing BCS instabilities, even at half- filling 48]. For sufficiently large t% there is a crossover to 
a regime where the BCS processes eventually dominate, signaling a superconducting instability. Furthermore, because 
the Fermi surface is imperfectly nested, the growth of both BCS and AF couplings weakens. Further increasing t% 
eventually destroys nesting of the Fermi surface altogether and both types of divergences are suppressed. Three cases 
illustrating the crossover are shown in Fig. 1151 

5.0 



4.0 
3.0 



2.0' 
1.0 

T).0 2.0 4.0 6.0 8.0 10.0 

scaling parameter / 

FIG. 15: At half-filling, spin-wave instability occurs for £2 = 0, but as ti increases, the BCS instability wins over. This is due 
to imperfect nesting. As ti is increased further, the nesting is destroyed and both divergences are suppressed. The hopping ii 
is chosen such that t\ + 1% = 1. 

As further increases in the diagonal hopping £2 suppress the d x 2_ y 2 BCS channel, this channel diminishes relative 
to other subdominant BCS channel with different symmetries, for example d^j,- or p-wave, as shown in Fig. EH These 
other channels, however, are all repulsive and hence do not lead to BCS instabilities by themselves. 

We note that the Hubbard model on the anisotropic triangular lattice has also been studied using the the random- 
phase approximation^^ and the fluctuation-exchange (FLEX) approximation!^ 0H1 E3 ■ The d-wa.ve superconduct- 
ing instability was found to be dominant for a large range of values of £2/(^1 + £2) interpolating between the square 
lattice and the isotropic triangular lattice. 
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FIG. 16: Dominant BCS channels for ti = 0.9 and t% = 0.1, at half-filling and at scaling parameter I — 5. These channels are 
both repulsive (with eigenvalues of 0.19525 and 0.162 respectively). The attractive d x 2_ y 2 BCS channel has a small eigenvalue 
of -0.268. 



E. Isotropic Triangular Lattice At Half-Filling 

For the special case of the isotropic triangular lattice (t\ = ti = 0.5) at half-filling, the weak-coupling RG flows do 
not show any BCS instabilities. The dominant BCS channels, d x 2_ y 2 1 d xy and p, are all repulsive. These channels are 
depicted in Fig. I17f A1. ITTT B) and !17f CL Fig. H7f D) shows the first attractive channel that develops, but the rapid 
oscillations in the effective potential, and its small size, indicate that the calculation is not reliable. 

In Fig. El the dominant AF and BCS channels are compared. Neither channel shows strong divergences, but the 
AF channel is significantly larger than the BCS channel. Thus there are signs of re-entrant antiferromagnetic long- 
range order. We speculate that there exist four different regions as the isotropy parameter t^j (t\ -M2) changes form 
(square lattice) to 1 (decoupled chains). These phases are shown in Fig. El At £2 = the system exhibits long-range 
antiferromagnetic Neel order (LRO) with ordering vector Q = (tt,tt). But this long-range order is suppressed by 
turning on £2- Instead d x 2_ y 2 BCS instabilities dominate and only short-range order (SRO) occurs. Both AF and 
BCS instabilities are suppressed as t-x is further increased and the nesting is completely eliminated. Nevertheless, the 
AF coupling remains significantly larger than the BCS coupling. In the strong-coupling limit U — > 00 the model can be 
mapped onto a nearest-neighbor spin-1/2 Heisenberg antiferromagnct which of course is insulating. On the isotropic 
triangular lattice this antiferromagnet exhibits long-range AF orderj5j||j^| with ordering wave vector Q — (An/ 3,0). 
The question of whether or not our weak-coupling analysis can describe this strong-coupling limit is tantamount to 
asking whether or not one or more intermediate-coupling fixed points intervene between the repulsive weak-coupling 
fixed point that is accessible in our RG analysis, and the attractive strong-coupling fixed point which describes the 
antiferromagnetic insulator. Finally as t% becomes larger than t\, the chains begin to decouple. In the extreme limit 
of independent chains there can be no long-range spin order, as the Mermin- Wagner theorem tells us that continuous 
symmetries cannot break in 1+1 dimensions. 

Curiously, four such regions, LRO / SRO / LRO / SRO, have also been identified in the corresponding strong- 
coupling Heisenberg antiferromagnet with two exchange couplings J% and J2. The phase diagram of this model 
has been studied via a straightforward 1/S expansion [Sfl IHfj. a series expansion method [57j. and a large- N Sp(N) 
solution [5J|. All three methods find two regions of long-range order: near the limit of a square lattice (J2 = 0) and 
near the isotropic point (J2 = J\). It is remarkable that our weak-coupling RG analysis shows similar behavior. 

The Hubbard model on the triangular lattice has been studied at intermediate values of the interaction strength 
within the Hartree-Fock approximation 59, 60] and within the slave-boson method [6Ul6al63| . A Mott-Hubbard metal- 
insulator transition is found to occur at a relatively large value of U c (U c = 5.27t for the Hartree-Fock calculation [5^ 
and U c — 7.23i or 7.68i from the slave-boson calculations [(HI EH). At a smaller value of U = U c \ there is also a 
continuous transition from a paramagnetic metallic phase to a metallic phase with incommensurate spiral order. The 
Hartree-Fock calculation yields U c \ — 3.97i and the slave-boson calculation gives U c \ = 6.68i for this transition. Signs 
of re-entrant AF order in our weak-coupling RG calculation, namely the relatively large size of the AF channels in 
comparison with the BCS channels, are broadly consistent with this picture, as AF tendencies can be a precursor to 
a transition to the insulating state. 
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FIG. 17: BCS channels for the isotropic triangular lattice ti = ti = 0.5, at scaling parameter I = 8.4. Channels with d x i_ y i 
(A), d xy (B) and p (C) symmetries all appear as repulsive channels (with eigenvalues 0.0378, 0.035 and 0.057, respectively). 
The largest attractive BCS channel is shown in (D), but it has a very small coefficient (-0.020) and the rapid oscillations suggest 
that the calculation is not accurate at this point. In the calculation, the Fermi surface was divided into M — 24 patches, instead 
of just 16, to improve the accuracy in the higher wave channels. 
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FIG. 18: Flow of dominant AF and BCS channels for the case of the isotropic triangular lattice. 
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FIG. 19: Four regions of the U > phase diagram, LRO / SRO / LRO / SRO, identified as t2/{t\ + £2) varies from (square 
lattice) to 1 (decoupled chains). 



IV. CONCLUSION 



Hubbard models have received extensive study in the context of high-T c superconductivity. We have reproduced 
the well-known result that there is an AF instability at half-filling on the square lattice with repulsive on-site Coulomb 
interaction and nearest-neighbor hopping. Furthermore, upon doping the system away from half filling, a crossover 
to a BCS regime with d x 2_ y 2 pairing symmetry occurs as expected. We have shown that it is important to retain 
subleading, formally irrelevant, corrections to the RG flows when the bare interaction is repulsive and the Fermi 
surface is nearly nested. We also studied another way of triggering a BCS instability. Keeping the system at half 
filling, but introducing the diagonal hopping £2 as shown in Fig. ^ along one of the two diagonals, breaks up perfect 
nesting. Corresponding magnetic frustration kills the spin density wave, and Cooper pairing dominates, at least if ti 
is not too large. This result suggests that superconductivity can occur in a model of strongly correlated electrons, 
even at half-filling. We emphasize that stripes are not expected to play a role here; even moderate on-site Coulomb 
repulsion should inhibit charge segregation at half-filling. 

The half-filled Hubbard model on the anisotropic triangular lattice has been proposed as a minimal description 
of the conducting layers of the k-(BEDT-TTF)2X organic superconductors. It is important to establish the pairing 
symmetry of the superconducting state of these materials. Our theoretical results predict pairing of the d x 2_ y 2 type 
and in fact signs of such order have been seen experimentally. We did not find any evidence of spontaneous time- 
reversal symmetry breaking. Pairing symmetry of the type d x 2_ y 2 ±id xy would occur if an attractive d xy channel arose 
in addition to the d x 2_ y 2 channel. We find that attractive d xy channels neither occur when next-nearest-neighbor 
hopping t' is included on the square lattice, nor when non-zero hopping ti is turned on along one of the two diagonal 
directions. 

Finally, we made contact with previous work on Hciscnbcrg antifcrromagnets on the anisotropic triangular lattice. 
Our weak-coupling RG calculation shows AF tendencies in two separate regimes - tendencies which seem to corre- 
spond with the two AF ordered phases found previously at large- U . In particular the portion of the phase diagram 
between the isotropic point (ii = £2) and decoupled chains (ti = 0) is the relevant region for the layered CS2CUCI4 an- 
tiferromagnet insulator material. The competition we found between antiferromagnetic order and spin-liquid behavior 
in our RG calculation may be consistent with the observed ease by which spin order is destroyed in the compound |27|. 
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